rm(list = ls())

library(tidyverse)

dat <- readRDS('data/figA1.rds')

## Plot 

p1 <- ggplot(dat, aes(x = n_outlets_now_unique, y = n)) + 
  theme_bw() + 
  geom_bar(stat = 'identity') + 
  geom_vline(xintercept = median(dat$n_outlets_now_unique, na.rm = T),
             linetype = 'dotdash') + 
  scale_x_continuous(breaks = seq(0, max(dat$n_outlets_now_unique), 1)) + 
  labs(x = 'Number of outlets in county',
       y = 'Number of observations\n(county-year)')

p1 

